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Search for unbound 15 Be states in the 3n+ 12 Be channel 
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Background: 15 Be is expected to have low-lying 3/2 + and 5/2 + states. A first search did not observe the 
3/2 + [A. Spyrou et al, Phys. Rev. C 84, 044309 (2011)], however, a resonance in 15 Be was populated in a 
second attempt and determined to be unbound with respect to 14 Be by 1.8(1) MeV with a tentative spin-parity 
assignment of 5/2 + [J. Snyder et al., Phys. Rev. C 88, 031303(R) (2013)]. 

Purpose: Search for the predicted 15 Be 3/2 + state in the three-neutron decay channel. 

Method: A two-proton removal reaction from a 55 MeV/u 4, C beam was used to populate neutron-unbound 
states in ls Bc. The two-, three-, and four-body decay energies of the 12 Bc + neutron(s) detected in coincidence 
were reconstructed using invariant mass spectroscopy. Monte Carlo simulations were performed to extract the 
resonance and decay properties from the observed spectra. 

Results: The low-energy regions of the decay energy spectra can be described with the first excited unbound 
state of 14 Be ( E x = 1.54 MeV, E r = 0.28 MeV). Including a state in ls Be that decays through the first excited 
14 Be state slightly improves the fit at higher energies though the cross section is small. 

Conclusions: A ls Be component is not needed to describe the data. If the 3/2 + state in 15 Be is populated, the 
decay by three-neutron emission through 14 Be is weak, < 11% up to 4 MeV. In the best fit, ls Be is unbound with 
respect to 12 Be by 1.4 MeV (unbound with respect to 14 Be by 2.66 MeV) with a strength of 7%. 

PACS numbers: 27.20+n,21.10.Dr 


The study of neutron-unbound states in atomic nuclei 
has been instrumental in probing the neutron drip-line 
[I]. Improvements in experimental techniques continue 
to push the discovery of nuclei farther from the valley 
of stability. Recent measurements of nuclei unbound by 
two neutrons include the first observations of 10 He [2j], 
13 Li (3|, 16 Be [i|, and evidence for 26 0 

The beryllium isotopes have been a fertile ground to 
study neutron-unbound states because of the availability 
of beams and low level densities. The heaviest beryllium 
isotope with a known bound state is 14 Be I6| (though 
13 Be is unbound) and the heaviest isotope measured is 
the two-neutron unbound 16 Be 0]. Between these two 
isotopes lies 15 Be, whose ground state still has yet to 
be confirmed. Shell model calculations have predicted 
low-lying 3/2 + and 5/2 + states located near each other 
[Jii, .1]. The first attempt to experimentally observe 15 Be 
was performed using a two-proton removal from a 17 C 
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beam j§}. The ground state of 17 C has been shown to 
have a 3/2 + spin and parity [9j and the reaction used 
was expected to remove two p-shell protons while leaving 
the neutrons in their initial configuration. Very few 14 Be 
fragments were observed in coincidence with the detected 
neutrons. Therefore it was concluded that the ground 
state of 15 Be is likely to be located at a higher energy 
than the first unbound excited state in 14 Be. This state 
at E r = 0.28 MeV would serve as an intermediate state 
for decay to 12 Be. A second experiment was performed 
using a neutron stripping reaction on a CD 2 target with 
a 14 Be beam to populate 15 Be 0- A 15 Be resonance 
was observed with a one-neutron decay energy of 1.8(1) 
MeV, a width of 575(200) keV and 1 = 2. The result 
of this experiment alone does not answer the question of 
the location and the properties of the 15 Be ground state, 
however, since the 3/2 + state could be above the first 
excited state in 14 Be and still be lower than the observed 
resonance in 15 Be. 

Shell model calculations predict a large spectroscopic 
overlap (1.27 for t = 2) between the predicted 3/2+ state 
and the first excited 2 + state in 14 Be Q. This calcula¬ 
tion, paired with the non-observation of a low-lying state 
in the two-body decay of 15 Be [§}, support the idea that 
the 3/2 + state could be observed in the three-neutron 
decay channel after decaying through the first neutron- 
unbound 2 + state in 14 Be. This is shown schematically 
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FIG. 1. (Color online) Select low-lying levels in beryllium 
isotopes from [lol - [T^ |. The dashed blue arrow represents the 
one-neutron decay from ls Be to 14 Be 0 and the solid red 
arrows represent the suggested decay path in § searched for 
in this work. The red-shaded box represents the range of 
energies that were simulated. 
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FIG. 2. (Color online) The detected fragments deflected by 
the Sweeper magnet are shown by their mass-to-charge ratio. 
The black histogram represents fragments detected in singles 
mode and the filled red histogram shows fragments detected 
with a neutron coincidence. Two vertical dashed lines show 
where 14 Be fragments should appear. This figure is adapted 
from Ref. |{J. 

in Fig. []]by the solid red lines. In the original analysis 
fU, few 14 Be fragments were observed while a large num¬ 
ber of 12 Be fragments were detected, as shown in Fig. 
[2j In the present work, we analyze the possible 3n+ 12 Be 
channel from the 2 p removal experiment in Ref. j§] by 
reconstructing decay energies with multiple neutrons in 
coincidence with those 12 Be fragments. 

The experiment was performed at the National Su¬ 
perconducting Cyclotron Laboratory at Michigan State 


University. A 55 MeV/u 4 'C beam was produced from 
a 22 Ne primary beam. This beam was focused onto 
a beryllium target where neutron-unbound states were 
populated. The emitted neutrons were detected with the 
Modular Neutron Array (MoNA) [T§]. The fragments 
were deflected by the Sweeper dipole magnet 14] into 
charged-particle detectors. The measured four-momenta 
of the neutron(s) and fragment allowed the reconstruc¬ 
tion of the decay energies by invariant mass spectroscopy. 
More experimental details can be found in Ref. Q. 

The multiplicity, or number of scintillator bars with 
a valid signal in a given event, is a useful parameter in 
searching for a state decaying by multiple neutron emis¬ 
sion. In scenarios where unbound states decay by one 
neutron, the multiplicity distribution is most likely to 
peak at one and decrease with each higher multiplicity 
value. This distribution is a result of scattering to mul¬ 
tiple bars. When resonances decay by more than one 
neutron there is a greater probability of detecting multi¬ 
ple neutrons. This shifts the distribution toward higher 
multiplicities. To search for decays of 15 Be to 12 Be+3n, 
decay energies were reconstructed using one, two, and 
three neutrons with a coincident 12 Be fragment in a time- 
ordered manner. In addition, restrictions were applied to 
the spectra to enhance or discriminate certain features. 
Applying a multiplicity-one restriction to a two-body de¬ 
cay energy spectrum enhances the signal from single neu¬ 
tron decays. To reduce cross-talk from multiple scatter¬ 
ing between scintillator bars in MoNA, causality require¬ 
ments were applied to the three-body decay energy spec¬ 
trum 0J, fim . The causality requirements used here 
restrict the first two hits to have velocities greater than 
10 cm/ns and the distance between hits to be greater 
than 30 cm. While many valid two-neutron events are 
removed, the ratio of true two-neutron events to multiple 
scattering events is greatly improved. There were insuf¬ 
ficient statistics to apply any causality cut to the four- 
body spectrum, thus it is dominated by cross-talk from 
one and two neutrons scattering in MoNA. Six spectra 
constructed from experimental data are shown in Fig. [3] 
the two-, three-, and four-body decay energy spectra, the 
multiplicity-one gated two-body spectra, the three-body 
spectra with causality cuts, and the multiplicity distri¬ 
bution. 

Monte Carlo simulations, which take into account ex¬ 
perimental acceptances, resolutions, and efficiencies were 
performed to provide an interpretation of the data jl8j . 
To describe the data presented in this work, previously 
observed states in 13 ’ 14 Be 0i were included in the 
simulations. These states are shown in Fig. [T] The en¬ 
ergies and widths were fixed, only the relative strengths 
were free parameters. In addition to the simulation of 
these known states, the simulation of a new state in 15 Be 
with freely varying energy, width, and relative strengths 
was included. The most sensitive parameter for this state 
was the decay energy which was varied from 50 keV to 4 
MeV. 

The simulated data are shown on top of the experi- 
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FIG. 3. (Color online) From left to right in the top row are the two-, three-, and four-body decay energies. From left to right 
in the bottom row are the multiplicity equal to one two-body decay energy, the three-body decay energy with causality cuts, 
and the multiplicity. The black crosses represent the experimental data. The dashed blue line represents the best fit if no 15 Be 
is included. In this fit the 0.28-MeV state in 14 Be is directly populated. The best fit is shown by the solid red line, which is 
the sum of the solid purple, 15 Be at 1.4 MeV decaying through the 14 Be 0.28-MeV state and dashed green lines made up of the 
higher lying 14 Be states. 


mental data in Fig. [3j The experimental decay energy 
spectra feature a low-energy peak. The energy and width 
of this peak, taking into account the experimental res¬ 
olution, can be described by the first excited neutron- 
unbound 2 + state in 14 Be. This state is unbound with 
respect to 12 Be by 0.28 MeV [lj|. The higher-energy 
components of the spectra can be reproduced by higher- 
lying states reported in Ref. [l2| (though, these param¬ 
eters are not unique solutions to the fit). The 0.28-MeV 
state accounts for a: 11% of the total strength up to 4 
MeV in the best fit when populated directly (no 15 Be 
component). The sum of the 0.28-, 2.28-, and 3.99-MeV 
states is shown by the dashed blue line in Fig. [3] 

A fit was also performed including a state in 15 Be as 
indicated by the solid red line in Fig. [3] This fit is made 
up of the 15 Be state decaying through the 0.28-MeV 14 Be 
state (solid purple line) and higher-lying states in 14 Be 
(dashed green line). The y 2 was calculated by comparing 
the simulated data to the experimental data in all six 
spectra. Using the y 2 value, the best fit (denoted by the 


star) and the one-, two-, and three-cr limits as functions of 
the strength of 15 Be and the 15 Be decay energy are shown 
in Fig. [IJ To obtain ratios below the 0.28-MeV state in 
14 Be, four-body decays without intermediate states were 
simulated. The strength of this state for the best fit 
is ss 7% and the ler limit is at ss 11%. The fit also 
includes 1.1% direct population of the 14 Be state at 0.28 
MeV. The simulations are insensitive to the width of the 
included 15 Be state. For the fit shown, an £ = 2 resonance 
with a total width of 500 keV (arbitrarily chosen due to 
insensitivity) was used. As the energy of the state is 
increased, the fit approaches that of the 0.28-MeV state 
in 14 Be being directly populated instead of acting as an 
intermediate state in the decay process. 

The overall cross section for populating 12 Be was cal¬ 
culated from the singles data (events not requiring a coin¬ 
cident neutron hit), shown in Fig. [2j to be 4.6 ±1.1 mb. 
The ratio of 15 Be to the total extracted from the sim¬ 
ulated decay energy spectra then implies that the two- 
proton removal cross section decaying to 12 Be is smaller 
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FIG. 4. (Color online) The one- (red), two-(blue), and three- 
a (green) limits as a function of the ratio of 15 Be to the total 
counts and the 15 Be decay energy above 12 Be. The best fit is 
indicated by the star. 


than 0.5 mb. This upper limit is a factor of two smaller 
than the calculated two-proton removal from 17 C of 0.99 
mb 0. In addition, it would also suggest that the cross 
section for populating 14 Be directly in a 2pln removal 
reaction would be ~ 4 mb, which would be larger than 
the two-proton removal cross section. Calculations pre¬ 
dict the cross section to be an order of magnitude lower 
for the 2pln removal process compared to the 2p removal 


To resolve this discrepancy, a different process con¬ 
tributing to the overall cross section was considered. Re¬ 
cently, Sharov et al. suggested that direct a removal 
might contribute to the population of 10 He in the reac¬ 
tion of 14 Be on 9 Be at 59 MeV/u [2lj]. The correspond¬ 
ing process in the present reaction would lead directly to 
13 Be which then would populate 12 Be by the emission of 
a single neutron. Similar evidence for such a process has 
also been observed in the analysis of the break-up data 
of 14 Be populating 10 He p2j . 

In order to test for this possibility, simulations were 
performed where two reaction mechanisms were included: 


two-proton removal to 15 Be and a removal to 13 Be. Ex¬ 
perimentally, 2p2n and a removals were indistinguish¬ 
able. However, the cross section for 2p2n removal is ex¬ 
pected to be even smaller than the 2 p removal. For the 
two-proton removal, a single resonance in 15 Be was as¬ 
sumed to decay via the 280-keV state in 14 Be and the 
a removal was allowed to populate any of the three pre¬ 
viously measured resonances in 13 Be [11], shown in Fig. 

[I] Peaks at similar energies have been measured in other 
previous works such as Refs. [23l - [25i ]. Free parameters 
were the energy and width of the state in 15 Be, as well as 
all the relative strengths. The resulting best fit is shown 
in Fig. [5j In this fit, the energy of the 15 Be state is 1.6 
MeV. Based on the strengths of the states, the cross sec¬ 
tion for a removal would then be ~ 4 mb. Such a large 
cross section favors a direct a removal over a two-proton 
removal indicating a small two-proton removal cross sec¬ 
tion. 


This work searched for the unobserved 3/2 + state in 
15 Be. Based on comparisons of simulated data to ex¬ 
perimental data, a 15 Be state is not required. A fit with 
2pln removal to directly populate unbound states in 14 Be 
is able to reproduce the decay energy and multiplicity 
spectra. If a state in 15 Be is populated, the state is un¬ 
bound to 12 Be by 1.4 MeV with a strength < 11% deter¬ 
mined by the minimum y 2 . Based on cross section cal¬ 
culations from the data for this scenario, the 2 p removal 
cross section would be nearly an order of magnitude less 
than the 2pln removal cross section. This disagrees with 
theoretical calculations which predict the 2pln removal 
cross section to be an order of magnitude less than the 2 p 
cross section. An alternative approach using simulated 
2 p removal to 15 Be which decay by three neutrons and 
a removal to 13 Be which decay by one neutron to 12 Be 
was also able to describe the data. In this solution the 
a removal is the dominant reaction and again the cross 
section for the population of 15 Be is small. 
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FIG. 5. (Color online) The six spectra are the same as Fig. 02 however, with a different fit. The solid black line represents the 
best fit made by summing the components of a state populated by two-proton knockout to 15 Be at 1.6 MeV decaying through 
the 0.28-MeV state in 14 Be (red dashed line) and a removal to 13 Be at 0.40 MeV (dotted-dashed blue line), 0.85 MeV (solid 
orange line), and 2.35 MeV (solid light blue line). 
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